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ABSTRACT 

We present atmospheric models of red giant stars of various metallicities, 
including extremely metal poor (XMP, \M\ < —3.5) models, with many chemical 
species, including, significantly, the first two ionization stages of Strontium (Sr ) 
and Barium (Ba ), treated in Non- Local Thermodynamic Equilibrium (NLTE) 
with various degrees of realism. We conclude that 1) for all lines that are useful Sr 
and Ba abundance diagnostics the magnitude and sense of the computed NLTE 
effect on the predicted line strength is metallicity dependent, 2) the indirect 
NLTE effect of overlap between Ba and Sr transitions and transitions of other 
species that are also treated in NLTE non-negligibly enhances NLTE abundance 
corrections for some lines, 3) the indirect NLTE effect of NLTE opacity of other 
species on the equilibrium structure of the atmospheric model is not significant, 
4) the computed NLTE line strengths differ negligibly if collisional b — b and 
b — f rates are an order of magnitude smaller or larger than those calculated with 
standard analytic formulae, and 5) the effect of NLTE upon the resonance line of 
Ba II at 4554.03 A is independent of whether that line is treated with hyperfine 
splitting. As a result, the derivation of abundances of Ba and Sr for metal-poor 
red giant stars with LTE modeling that are in the literature should be treated 
with caution. 



Subject headings: stars: atmospheres, late-type — stars: abundances — line: for- 
mation 
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1. Introduction 

Generally, the time development of the neutron capture (n-capture) element (Z > 30) 
abundances over the galaxy's history encodes information about the astrophysical sites of 
the various n-capture processes, such as the "classical" rapid process (r-process), and the 
weak and main slow processes (s-pro cesses), and their relation to the sites of light metal 
(Z < 30) nucleosynthesis (see Cowan & Thielemann (2004) for a general review). High 
quality spectroscopy of metal-poor and extremely metal-poor (XMP) red giants in the halo 
component of the Milky Way galaxy has allowed the determination of n-capture element 
abundances at very early times in the formation of the galaxy (see James et al. (2004) 
(hereafter: JFBCGS04), for a recent example). In particular, the abundances of Sr and 
Ba are significant because, if their formation is understood correctly, they may indicate the 
relative rate of nucleosynthesis of light and heavy n-capture elements, and, by extension, 
the relative importance of the weak and main s-processes in the chemical evolution of the 
galaxy (see JFBCGS04 and references therein). 

Almost all large scale abundance surveys of old red giants to date have been carried out 
with atmospheric models and synthetic spectra computed under the assumption of Local 
Thermodynamic Equilibrium (LTE), either with MARCS models (Edvardsson et al. 1993) 
or ATLAS9 models (Kurucz 1994). Generally, previous investigation has indicated that the 
equilibrium state, and, thus, the line opacity, of many elements may deviate from LTE, 
particularly in giant atmospheres where gas densities, and, hence, thermalizing collision 
rates, are low (see Short & Hauschildt (2003) for example). Furthermore, it has been found 
that the extent of NLTE effects depends on the overall metallicity, [^] of the atmosphere. 
Mashonkina et al. (1999), (hereafter: MGB99), carried out Non-LTE (NLTE) modeling of 
the Ba II lines for main sequence and turn-off stars using a complete linearization method to 
solve the NLTE statistical equilibrium (SE) rate equations given an LTE atmospheric model 
structure. The study of MGB99 is an important pioneering calculation, and concluded 
that NLTE effects affected the inferred abundance of Ba by as much as 0.2 dex for stars 
of [^] < —2. However, there are three things about the study of MGB99 that we draw 
attention to: 1) it is not self-consistent in that the model structure adopted does not take 
into account NLTE effects on the many other species, such as Fe, whose opacities determines 
the radiative equilibrium structure, 2) it is not self-consistent in that other species that may 
have overlapping transitions with Ba are treated in LTE, 3) because of the inefficiency of the 
complete linearization method the Ba II atomic model is incomplete, having 35 levels and 
(judging by the Grotrian diagram in their Fig. 1) 28 bound-bound (b — b) transitions. It is 
not clear that these limitations affect the results significantly without performing numerical 
experiments with the same models. However, as described in Section 2, our NLTE modeling 
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should be more realistic on all three of these points. 

The atmospheric and spectrum modeling with the PHOENIX code for the mildly metal- 
poor red giant Arcturus (a Boo, HD124897, spectral type K2III) presented by Short & 
Hauschildt (2003) included many thousands of lines of many species in NLTE, but did 
not include Sr or Ba . We have recently added to the NLTE atmospheric and spectrum 
synthesis code PHOENIX the ability to compute both of these species in NLTE. The purpose 
of the present investigation is to determine the effect on the computed Sr and Ba lines, 
and, hence, on the inferred abundances, of three types of NLTE consideration: 1) direct 
NLTE effects on the Sr and Ba equilibrium and line formation, 2) indirect NLTE effects on 
the atmospheric structure from NLTE effects on the opacity of many other chemical species, 
and 3) the indirect NLTE effect on the computed Sr or Ba state by overlap of Sr and Ba 
transitions with transitions of other species also in NLTE. Moreover, we investigate how the 
NLTE effects vary with overall value of [A], including XMP values. Here, the investigation is 
carried out for the most readily observable metal poor stars in the galactic halo: red giants. 
In Section 2 we describe the model calculations, in Section 3 we present our results, and we 
re-iterate our main conclusions in Section 4. 



2. Modeling 

2.1. Stellar parameters 

To represent a typical XMP red giant of the type that has been observed in the galactic 
halo we adopt the following stellar parameters: a T eff value of 4800 K, a logg value of 
1.5, a value for the microturbulent velocity dispersion, £t, of 2.0 km s -1 , and a convective 
mixing length parameter of 2 pressure scale heights. The models are computed with spherical 
geometry with a radius at the ri 2 ooo surface of 1.6 x 10 12 cm ( 25 R©), where T12000 is the 
optical depth due to continuous opacity at a standard reference wavelength of 12000 A. 
These parameters were chosen to be approximately representative of two XMP red giants 
recently described by Cayrel et al. (2004): CD -38:245 and CS 22949-037. 

Abundances: A set of models of these parameters was generated with a scaled solar 
abundance distribution spanning a range in |A] value from -1 (mildly metal poor) to -5 
(XMP), with an interval of -1. We adopt the solar abundances of Grevesse et al. (1992). 
In all cases the values of [^] and [^] equal the value of [A]. The use of a scaled solar 
abundance distribution is known to be unrealistic for Pop II, and especially extreme Pop II, 
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stars. However, there is a range of abundance distributions that has been found for XMP 
stars, so no one choice would represent all of them. Furthermore, the purpose of our study 
is to investigate differential effects of [g] value on the NLTE spectrum formation of Sr and 
Ba , and, to begin with, we suppose that the abundance distribution will be unimportant to 
first order for such differential effects. 



2.2. PHOENIX 

PHOENIX makes use of a fast and accurate Operator Splitting/Accelerated Lambda It- 
eration (OS/ALI) scheme to solve self-consistently the radiative transfer equation and the 
NLTE statistical equilibrium (SE) rate equations for many species and overlapping transi- 
tions (Hauschildt & Baron 1999). In our fully self-consistent models (NLTEp u ii models, de- 
scribed below) we include a total of 39 species in the set of multi-level NLTE SE equations, 
including the lowest one to three ionization stages of most light metals up to and including 
the Fe group, as well as the first two ionization stages of Sr and Ba . We have only in- 
cluded ionization stages that are significantly populated at some depth in the atmosphere 
of a late- type star, as computed in radiative-convective equilibrium (ie. no chromosphere or 
transition region). Therefore, no ionization stages higher than III are treated in NLTE for 
any species. Short & Hauschildt (2005) contains details of the number of ionization stages, 
E levels, and b — b and bound-free (b — f) transitions included in the NLTE PHOENIX SE 
calculation for late-type stars, as well as the sources of atomic data and the formulae for 
various atomic processes. 

2.3. Degrees of realism 

Table 1 shows the four degrees of realism in our modeling: 

1) NLTEpuu: the most self-consistent NLTE models with 35 other species in addition to 
Sr and Ba I and II included in NLTE in the calculation of both the atmospheric T 
and N e structure and the equilibrium state of Sr and Ba . Included in NLTE SE are 
species such as Fe I and II that are known to have an important effect on the computed 
radiative equilibrium T structure. These models account for the indirect effects of both 
the NLTE opacity of many species on the atmospheric T and N c structure, and multi-level 
NLTE effects caused by transitions of other species in NLTE that overlap with Sr and Ba 
transitions, as well as the direct effects of NLTE on Sr and Ba themselves. 
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2) NLTE Sr +B a : only Sr and Ba treated in NLTE SE with the NLTE FuU T structure of 1). 
These models account for the effects of NLTE opacity of other species on the atmospheric 
structure, and for NLTE effects on Sr and Ba , but not for the effects of NLTE on transitions 
of other species that overlap those of Sr and Ba . 

3) NLTElte: deliberately inconsistent models in which the state of Sr and Ba is computed 
in LTE (ie. no species in the NLTE SE), but with the NLTEp u ii atmospheric structure of 1). 
The purpose of these models is to distinguish between effects on the Sr and Ba lines that 
are caused by direct NLTE effects on the Sr and Ba equilibrium and those that are caused 
by the NLTE atmospheric structure. 

4) LTE: fully LTE models with the atmospheric T and N c structure and the state of all 
species, including Sr and Ba , computed in LTE. These are the simplest models and are 
equivalent to the MARCS and ATLAS9 models that are largely used in abundance analyses. 

2.4. Sr and Ba I and II 

Table 2 contains details of the numbers of E levels and b — b transitions of Sr and 
Ba I and II included in the NLTE SE rate equations. All levels not designated as 
"Sum" are fine structure (J) levels, and were included separately in the NLTE SE rate 
equations. Energy level values, log gf values, and damping parameters have been taken 
from Kurucz (1994). Because the purpose of our study is to investigate the differential effect 
of modeling treatments rather than to fit observed spectra, we have not performed a quality 
control investigation or fine-tuned the atomic data, except for the collisional cross-section 
perturbation analysis that follows. As a point of information, we note that many of the 
fine structure level designations for these species in the atomic data base of Kurucz (1994) 
are erroneous. Among the the levels included in our atomic models, all levels of Ba and 
select levels of Sr were fixed to agree with the designations given in the National Institute 
of Standards and Technology (NIST) Atomic Spectra Database (for Ba ) or Handbook of 
Basic Atomic Spectroscopic Data (for Sr ). Unfortunately, we were unable to find data for 
many of the Sr levels in any of the readily accessible databases, so were unable to check 
the designations for all the levels. As with all the species treated in NLTE, only levels of 
Sr and Ba connected by transitions of log gf value greater than -3 (designated primary 
transitions) are included directly in the SE rate equations. All other transitions of that 
species (designated secondary transitions) are calculated with occupation numbers set equal 
to the Boltzmann distribution value with excitation temperature equal to the local kinetic 
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temperature, multiplied by the ground state NLTE departure coefficient for the next higher 
ionization stage. 

3. Results and Discussion 

3.1. NLTEpuu vs LTE 

Fig. 1 shows the computed atmospheric T and N c structures for the NLTE Fu n and 
LTE models for all values of [^]. The NLTEp u ii models include much of the Fe group 
(the first two ionization stages of Ti , Mn , Fe , Co , and Ni ) in NLTE. These species, 
especially Fe , contribute much of the line opacity to the blue and UV bands and have a 
large effect on the computed radiative equilibrium of the atmosphere. Treating them in NLTE 
generally decreases T near the bottom of the atmosphere and increases it near the top of 
the atmosphere, thereby exhibiting a reduction in the classical "back-warming" effect of line 
opacity. However, the difference in NLTE and LTE T structures is negligible throughout 
much of the line forming region near logri2ooo ~ 1- See Short & Hauschildt (2005) and 
Anderson (1989) for a detailed discussion of thermal equilibrium in late- type atmospheres 
with massive scale NLTE. There is a rapid increase in the boundary temperature at the top 
of the model between [^] values of -1 and -3, and the model of [g] = —2 shows complicated 
behavior between logri 2 ooo values of -4 and -6. However, the [g] = —2 model has converged 
at all depths to a stable equilibrium, and its behavior is caused by [^] = —2 being a critical 
value for the thermal equilibrium in the outer atmosphere. 

Figs. 2 and 3 show computed line profiles for the lines used as abundance diagnostics by 
JFBCGS04, with all wavelengths in air: Ba II A 4554.03, 5853.69, 6141.73, and 6496.91, and 
Sr II A 4077.71 and 4215.52 computed with the most and least realistic modeling, NLTE Full 
and LTE, respectively. Selected atomic data from Kurucz (1994) for these lines are presented 
in Table 3, with the designations of the atomic fine structure levels fixed to agree with the 
NIST Atomic Spectra Database or the Handbook of Basic Atomic Spectroscopic Data. For 
each line we show the profiles for those models of [g] value for which that particular line is of 
detectable strength, yet on the linear part of the curve of growth (COG) (ie. unsaturated), 
or barely saturated, and thus maximally useful as an abundance diagnostic. These are also 
the lines for which NLTE effects, and the realism of their treatment, will make the greatest 
difference. Table 4 shows equivalent widths, W\, for these lines, computed with the various 
degrees of modeling realism. Again, for each line we only show W\ values for models of [g] 
value such that the line is detectably strong, yet unsaturated. 
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For both the Ba and Sr lines, the effect of NLTE as compared to LTE is to increase the 
predicted strength, W\, of the line for models of higher [g] value and decreases it for models 
of lower [g] value. We find that the critical [g] value varies from line to line throughout 
almost the full range of [g] values spanned by our models. Furthermore, for each line studied, 
the NLTE effect is significant at some value of [g]. For unsaturated lines, NLTE corrections 
to and [^] are proportional to log W\ in the sense that a model that predicts too large 
a value of log W\ because of inadequate physical treatment will cause the inferred value of 
[ff] or [if] to be proportionally too small, and visa versa. We note that these results are 
qualitatively the same as those of MGB99 in their study of main sequence and turn-off stars. 

Figs. 4 and 5 show for the NLTE Full models of [g] value equal to -1 and -4 the logarithm of 
the NLTE i?-level departure co-efficients, h — for each i?-level, i, of Sr and Ba stages I 
and II, and the ground state of stage III, where ri\ is the actual NLTE occupation number 
of level i and n* is the occupation number of the level as computed with LTE equilibrium 
relations for all levels and stages with the local value of T^ n and the LTE N e . For the 
models of [g] = —1, the b\ values for the ground level and the levels of lowest excitation 
energy of stage II remain close to unity throughout the atmosphere. This is expected 
because stage II is the dominant ionization state and so its total number density, most 
of which will be in the ground state, should be negligibly affected by departures from the 
LTE ionization equilibrium. However, for the model of [g] = —4, the ground state b\ value 
for Ba II drops below unity for T12000 < —4, and approaches 0.1 at the surface. The reason 
for the difference in the ground state b\ behavior in the models of [^] = — 1 as compared to 
those of —4 is that in the latter Ba II is no longer as dominant in the ionization equilibrium, 
but is rivaled by Ba III. For both [^] values, many levels of higher excitation energy in 
stage II are overpopulated with respect to LTE (b- x > 1) in the upper atmosphere due to 
the influence of collisional coupling to the ground level of stage III. The ground level of 
ionization stage III is overpopulated with respect to LTE in the outer atmosphere as a 
result of NLTE overionization of stage II (note that the ground level is the only level of 
stage III included in our SE equations). The ground state of stage I is underpopulated 
with respect to LTE in the upper atmosphere, again, as a result of NLTE over-ionization. 
In Fig. 8 we show the LTE and NLTE partial pressures of the three ionization stages of Ba 
for models of [jj] value equal to -1 and -4. We note that the results for the 6; values are 
consistent with the partial pressure results. 

From a comparison of Figs. 4 and 5 to Fig. 2 of MGB99 we note that for the models 
of [g] = —1 the results for the ground levels of Ba II and III are qualitatively similar 
to theirs for the Sun, where their ground state value of b\ for Ba II remains close to unity 
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throughout the atmosphere. The results for the ground level of Ba II for our [g] = —4 
model are also qualitatively similar to those of MGB99 for their model of G84-29, which has 
[g] = —2.6. The latter also shows a drop below unity in the upper atmosphere. However, 
these comparisons are not expected to be exact because out models correspond to giant stars, 
whereas their correspond to main sequence stars. 

3.2. NLTEfuii vs NLTEs r+ Ba, NLTElte and LTE 

Figs. 6 to 7 demonstrate the sensitivity of the results to the degrees of realism and 
completeness in the treatment of the NLTE problem by comparing the NLTEs r +Ba and 
NLTElte results to those of the NLTEpuii and LTE results. Note that we only carried out 
NLTEs r +Ba and NLTElte modeling for values of [g] equal to -1, -4, and -5, and in Figs. 6 
to 7 we only show profiles for those [jj] values for which that particular line is unsaturated. 

We find that in some cases (Ba II A 5853.69, 6141.73 and 6496.91 for [|] = -1) the 
less realistic NLTEs r +Ba models yield line profiles that differ significantly from those of the 
NLTEpuii model and are closer to the LTE profiles than the NLTEpuii profiles are. This 
indicates that these lines are indirectly affected by a significant amount by the non-local 
effect of overlapping transitions with other species that are also affected significantly by 
NLTE. In the case of these transitions of these species, the result of this "multi-species non- 
locality" is to enhance the deviation from the LTE results. Therefore, spectrum synthesis 
calculations in which only the species of immediate interest are treated in NLTE (Sr and 
Ba ) will underestimate NLTE corrections. This result points to the importance of a NLTE 
treatment for all species with significant line opacity for the correct NLTE solution for any 
one species. Particularly in late type stars where the spectrum is blanketed by 10 7 lines, 
there may be any number of other species with transitions that overlap one or more of those 
of Sr or Ba whose equilibrium state thus determines the supply of photons that are available 
for absorption by the Sr or Ba transition that is so overlapped. (Moreover, the equilibrium 
state of those other species may in turn be affected by yet other species that deviate from 
LTE whose transitions overlap those of those species!) The significance of this result bears 
careful consideration given the number of NLTE abundance determinations in the literature 
that are based on LTE "background" line opacity for all species other than the one under 
investigation. (Indeed, a way of understanding this result is that from the perspective of 
NLTE SE, there is no such thing as "background" opacity; the set of NLTE SE equations 
non- locally couples the states of all species.) 
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However, we caution that the magnitude of this multi-species non-locality effect for a 
particular line, or particular species, depends sensitively on the degree of wavelength overlap 
between lines of the species of interest that are important for the excitation equilibrium, and 
strong lines of other species treated in NLTE. The degree of overlap and the magnitude of 
its influence, in turn, depends on the accuracy of the atomic E level and oscillator strength 
(gf) data that is used to calculate the wavelengths and strengths, respectively, of the lines 
of overlapping species. The atomic data of Kurucz (1994) has been produced originally for 
the purpose of computing the broad-band solar flux level in LTE, for which only statistically 
correct line opacity is important, and is known to contain erroneous level energies and gf 
values. Therefore, the degree of the effect predicted here may be subject to revision, either 
upward or downward, if more accurate atomic data were used for all species in the NLTE 
calculation. Due to the non-locality of this effect, an exploration of the sensitivity of the 
effect to the perturbation of various E levels of various species is beyond the scope of the 
present investigation. 

In all cases the profiles computed with NLTElte modeling differ negligibly from those of 
LTE modeling. This indicates that the indirect effect of the NLTE opacity of other species 
on the atmospheric structure is much less important than the direct effect of NLTE SE on 
the Ba and Sr equilibrium and the less direct effect of overlapping transitions with other 
NLTE species. This result is consistent with the minor effect on the T structure of NLTE 
seen in Fig. 1. 

3.3. The role of collisional rates 

The greater realism of NLTE modeling comes at the expense of reliance on a greater 
number of physical input data. Among those that are crucial to determining the NLTE SE 
equilibrium and that are among the least well determined are the collisional cross-sections, 
cry, for b — b and b — f processes. The first E level above the ground state of Ba II has an 
excitation energy of 0.604 ergs, which corresponds to a T of 4700 K for an ideal monatomic 
gas. Therefore, collisional b — b transitions may be important for the excitation equilibrium 
of Ba II. In our modeling we use approximate values for e~ collisions calculated with the 
formula reproduced in Lang (1999) for forbidden transitions for all b — b transitions and of 
Drawin (1961) for b — f transitions. Note that, for consistency of treatment among b — b 
rates, and to reduce reliance of uncertain gf values, we do not use separate formulae for pairs 
of levels connected by permitted and forbidden transitions. Adding to the uncertainty in 
collisional rates is the possible role of collision partners other than electrons; MGB99 found 
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that collisions with H atoms play a detectable role in the NLTE equilibrium state of Ba II 
in main sequence and turn-off stars. 

To accommodate both of these considerations we perform a perturbation analysis by vary- 
ing all the b — b and b — f cry values for Ba II by a factor of 0.1 and 10. According to Fig. 
8, Ba I is at most 10" 3 times as abundant as Ba II and, thus, Ba I rates are not expected to 
affect the Ba II equilibrium. Moreover, there are no Ba I lines that are useful for abundance 
analysis. Therefore, we leave the Ba I rates unperturbed and vary separately both the b — b 
and b — f cry values for Ba II by a factor of 0.1 and 10. This approach assumes that any errors 
in cry values for either b — b or b — / processes are systematic among all transitions rather 
than random. However, a transition-by-transition perturbation analysis, even if limited to 
critical transitions, is beyond the scope of this investigation. The line opacity of Ba is not 
expected to affect the equilibrium atmospheric structure, so, with the cry values at their 
perturbed values, we re-converge the the Ba NLTE SE and recompute the spectrum using 
the NLTEpuu model structure. 

Fig. 9 shows the computed Ba II lines for the NLTE Fu ii model with various values for 
the cry collisional cross-sections. We find that perturbation of both the b — b and the b — f 
rates by a factor of 100 produces both departure coefficients and line profiles that differ 
negligibly from the those calculated with the standard approximations. This indicates that 
the excitation and ionization equilibria populations of Ba II is robust. In the case of the 
ionization equilibrium, the result is expected given that Ba II is the dominantly populated 
ionization stage. 

3.4. Hyperfine splitting of the Ba A4554.03 line 

The importance of hyper-fine splitting (HFS) of the Ba lines, especially the A4554.03 
line, has been well documented by McWilliam (1998), McWilliam et al. (1995), McWilliam 
& Rich (1994), and references therein, and these authors to various degrees present the 
formulae and data necessary for calculating the HFS components. HFS has the effect of 
de-saturating lines and thus significantly alters the inferred abundance from lines that are 
barely saturated and not yet heavily damped. For lines on the linear part of the COG, 
HFS changes the line profile, but not the value of W\. Because the hyper-fine sub-levels 
of a splitted level are relatively close to each other in E, they are expected to be strongly 
coupled by collisional transitions, despite such transitions being forbidden. Therefore, all 
such sub-levels should have the same b- x value as the unsplitted level, and we thus expect 
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that the effect of HFS on the profile of a line computed in NLTE will be the same as that 
for a line computed in LTE. However, for the sake of completeness, and because the HFS 
line profile is of interest, we investigate the effect of HFS in NLTE. 

We have modified our Ba II atomic model to account for HFS in the A4554.03 resonance 
line where the effect is known to be especially important. The structure of the PHOENIX 
code is such that it is more straightforward to represent the splitted sub-levels directly in 
the atomic model (the way nature does it) rather than by mimicing the effect of splitting 
by modifying the line profile function, 0„, as is suggested by, for example, Gray (1992). We 
deduced the E level splitting and the component log gf values from the component A and 
relative gf values given in Table 1 of McWilliam (1998). In stars, Ba may be significantly 
abundant over the range of atomic masses, A, from 134 to 138. The component gf values 
are weighted by the relative abundance of the isotopes of odd A value because only they 
have HFS. Because we are modeling metal-poor red giants, we follow McWilliam (1998) and 
Sneden et al. (1996), and adopt a pure r-process isotopic mix of 0:40:0:32:28 for isotopes of 
A value of 134, 135, 136, 137, and 138, respectively, which amounts to an odd to even mix 
of 72:28. 

Fig. 10 shows the A4554.03 line computed with the NLTEp u ii and LTE models with and 
without HFS in the treatment. We investigate the model of [^] value equal to -4 for which the 
line is clearly on the linear part of the COG. As expected, the effect of HFS is to broaden the 
line while preserving its strength. We note that HFS also makes the line profile asymmetric 
as a result of the asymmetry of the HFS components. Also as expected, the effect of NLTE 
on the profile computed with HFS is the same as that for the profile computed without HFS; 
the extent of the NLTE weakening of the line is independent of the E level splitting. 

4. Conclusions 

From comparison of LTE models with the most realistic NLTEp u ii models, we find that 
for all lines that are useful Sr and Ba abundance diagnostics the magnitude and sense 
of the computed NLTE effect on the predicted line strength is metallicity dependent, with 
stars of [jj] = — 1 to —2 showing NLTE strengthening of the lines and those of [^] = —4 to 
—5 showing NLTE weakening. These results are qualitatively similar to those of MGB99 for 
main sequence and turn-off stars and point to the need to account for NLTE effects when 
deriving spectroscopic Sr andBa abundances. 
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Ba : From Table 4, we see that among the most affected lines is the Ba A4554.03 
resonance ) for [|] = -4, for which log W A (NLTE Full )-log W A (LTE)=-0.14 (neglecting HFS), 
leading to a NLTE abundance correction of €nlte ~ ^lte — +0.14 (because the line is on 
the linear part of the COG). For Ba A5853.69 in models of [|] = -1, log T¥ A (NLTE Full )- 
logW A (LTE)=+0.04, leading to e NL TE - £lte = -0.04. 

Sr : Similar effects, though less pronounced, are found for the two Sr lines. We note that 
these lines are so strong that they are only on the linear part of the COG for the most metal 
poor of our models ([^] = —4 to —5). 

A less well studied NLTE effect is multi-species non-locality that arises from overlapping 
transitions that link the equilibrium state of one species with those of other species whose 
equilibrium also departs from LTE. By including most light metals up to and including 
much of the Fe group in NLTE (NLTEp u n models), we find that for some lines such overlap 
enhances NLTE abundance corrections non-negligibly, and should, in principle, be accounted 
for for accurate abundance determination. From Table 4 the largest such effect is for Ba 
A5853.69 for [|] = -1 for which log W^NLTEWHog W A (Sr+Ba)=+0.03. However, we 
reiterate that the magnitude of this effect is dependent on the accuracy of the atomic data 
of the overlapping NLTE species. Given the notorious innacuracy of such data, caution is 
required on this point. Another indirect NLTE effect is that of the opacity of other species 
that depart from LTE that are important for radiative equilibrium of the atmosphere. We 
find that the NLTE T structure is close enough to that of LTE in the line forming region that 
the effect on the computed Sr and Ba line profiles is negligible. We conclude that NLTE 
spectrum formation calculations with LTE atmospheric models will suffice for red giants. 
Freedom from the necessity to calculate large grids of NLTE atmospheric models for large 
NLTE abundance surveys will expedite such work. 

We find that the computed NLTE line strengths differ negligibly if collisional b — b and 
b — f rates are an order of magnitude larger or smaller than those calculated with standard 
analytic formulae. Because collisional cross-sections are notoriously uncertain, the robustness 
against perturbation in this case provides some assurance that the results are not artifact of 
erroneous atomic data. 

For completeness, we have also studies the effect of NLTE on the Ba A4554.03 resonance 
lines when HFS is taken into account. HFS for Ba has been extensively studied elsewhere 
and we reproduce the well-known result that HFS substantially widens the line and increases 
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the residual flux at line center. Of note here is that, as expected from the close collisional 
coupling among the splitted sub-levels, the relative effect of NLTE is the same with and 
without HFS. 
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Fig. 1. — Computed T and iV e structures of the atmospheric models as a function of the 
continuum optical depth at 12000 A. NLTE Fu ii: solid line; LTE: dashed line. 
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Fig. 2. — Selected computed Ba II lines for models of [g] values as indicated on each panel. 
Full NLTE model: solid line; LTE model: dashed line. Upper left: A4554.03 (resonance); 
upper right: A5853.69; lower left: A6141.73; and lower right A6496.91. Note that for each 
line we only show profiles for models of [g] value such that the line is of detectable strength, 
yet unsaturated or only barely saturated. 



-17- 




Fig. 3.— As for Fig. 2, but for Sr II. Left panel: A 4077.71; right panel: A 4215.52. 
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Fig. 4. — Departure coefficients of the NLTE statistical equilibrium solution for the model of 
[g] = —1. Left panels: Sr ; right panels: Ba ; upper panels: ionization stage I; lower panels: 
ionization stages II and III. The ground state coefficient is shown with a thick black line. 
The lighter the color of the line the higher the energy, E, of the level with respect to the 
ground state. The dashed line in the lower panels is the ground state departure co-efficient 
of stage III. 
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Fig. 5. — As for Fig. 4 except for the model of [g] = —4. 
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Fig. 6. — As for Fig. 2 with the addition of the lines computed with the NLTE Sr +Ba (dot- 
dashed line) and NLTE LT e modeling (crosses). 
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Fig. 7. — As for Fig. 3 with the addition of the lines computed with the NLTE Sr +Ba (dot- 
dashed line) and NLTE LT e modeling (crosses). 
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Fig. 8. — Partial pressures of Ba I (solid line), Ba II (dotted line), and Ba III (dashed line) 
for the LTE (dark line) and NLTE (light line) models. Models of [|] = -1 (upper panel) 
and [g] = —4 (lower panel). 
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Fig. 9. — Selected Ba II lines computed with the NLTEp u ii model with the rates from 
standard analytic approximations (see text) (solid line), with b — b rates decreased (dotted 
line) and increased (dashed line) by a factor of 10, and with b — f rates decreased (dot-dashed 
line) and increased (dot-dot-dot-dashed line) by a factor of 10. Note that the perturbation 
results are barely visible because they differ negligibly from the result with standard formulae. 
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Fig. 10. — As for Fig. 2 for the model of [g] = —4 with the addition of the model with 
hyper-fine splitting (HFS) (gray lines). 
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Table 1: Degrees of modeling realism. 

Model Species in Atmospheric 

NLTE SE structure 

NLTEFull Full set NLTE 

NLTESr + Ba Sr , Ba NLTE 

NLTELTE None NLTE 

LTE None LTE 
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Table 2: 


Sr and Ba 


I and II atomic models. 




Species 


No. £-leve 


Is No. b — b transitions 


Xi (eV) 


Sri 


52 


74 


5.70 


SrII 


32 


90 


11.03 


Bal 


76 


114 


5.21 


Ball 


51 


121 


10.00 
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Table 3: Selected atomic data for Sr II and Ba II lines of diagnostic utility. Data are from 
Kurucz (1994), except for level designations, which are from NIST. 



A A (Air) 


4i (s- 1 ) 


Xi (cm x ) 


Xi 


Lower level 


Upper level 


Ba II 


4554.033 


l.lle+08 


0.000 


21952.404 


[Xe }6s 2 S 1/2 


[Xe }6p 2 P° /2 


5853.675 


6.00e+06 


4873.850 


21952.404 


[Xe }5d 2 D 3/2 


[Xe }6p 2 P° /2 


6141.713 


4.12e+07 


5674.824 


21952.404 


[Xe }5d 2 D 5/2 


[Xe }6p 2 P° /2 


6496.898 


3.10e+07 


4873.852 


20261.561 


[Xe }5d 2 D 3/2 


[Xe WPli 2 


Sr II 


4077.71 


1.42e+08 


0.000 


24516.650 


[Kr }5s 2 S 1/2 


[Kr ]5p 2 P° /2 


4215.52 


1.27e+08 


0.000 


23715.190 


[Kt}5s 2 S 1/2 


[Kr ]5p 2 P° /2 
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Table 4: Equivalent widths, W\, for lines used in the abundance analysis of JFBCGS04 
under various modeling assumptions, shown for select [^] values. For unsaturated lines, the 
inferred abundance in [g] units is proportional to log 



H a (niA) 

Line [|] NLTEpuii NLTE Sr+B a NLTE LTE LTE 



Ball A4554.03 


-4 


48.4 


49.2 


66.6 


66.3 


Ball A5853.69 


-1 


132.2 


124.3 


121.8 


120.9 


Ball A6141.73 


-1 


182.7 


177.5 


172.7 


170.9 


Ball A6141.73 


-4 


16.1 


17.5 


17.6 


17.4 


Ball A6496.91 


-1 


179.3 


171.3 


167.7 


165.9 


Ball A6496.91 


-4 


11.8 


12.7 


12.8 


12.6 


SrII A4077.71 


-5 


47.2 


49.0 


55.0 


55.0 


SrII A4215.52 


-5 


32.6 


34.1 


39.0 


39.0 



